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^ Abstract: Current limit on the dark matter relic abundance may suggest that \n\ 

should be smaller than prediction in the minimal supergravity scenario (mSUGRA) 
for moderate m and m^j. The electroweak-ino parameter Mi,M 2 and are then 
much closer to each other. This can be realized naturally in the non-universal Higgs 
mass model (NUHM). Since the heaviest neutralino (xl) and chargino (xf) have sig- 
nificant gaugino components, they may appear frequently in the left-handed squark 
decay and then be detectable at the LHC In such a case, we showed that the hi- 
erarchy of Mi,M 2 and \fi\ can be determined. In the light slepton mass scenario 
with non- vanishing lepton-flavor violation (LFV) in the right-handed sector, NUHM 
with small corresponds to region of parameter space where strong cancellation 
among leading contributions to Br(fi — > c-f) can occur. We showed that determina- 
tion of electroweak-ino hierarchy plays a crucial role in resolving cancellation point 
of Br(n — > ej) and determination of LFV parameters. We also discussed test of the 
universality of the slepton masses at the LHC and the implications to SUSY flavor 
models. 
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1. Introduction 

Even though the Standard Model (SM) describes the Nature very well, there are a 
number of theoretical and phenomenological issues that the SM cannot give a plau- 
sible explanation. The notorious one is the hierarchy problem. In the absence of 
a symmetry to protect the mass of the Higgs boson from the radiative corrections, 
its natural value should be of order of UV cutoff scale. However, from the preci- 
sion electroweak measurements, the SM Higgs mass is lower than 193 GeV at the 
95% confidence level PJ. On the cosmo logical side, while the evidence of the dark 
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matter (DM) in the universe is now available from a variety of observational data 
000? the SM cannot provide a viable candidate for it. Therefore, it is 
convincing that the SM must be viewed as a low energy effective theory of the more 
fundamental theory. 

Supersymmetry (SUSY) is one among the most promising candidates of the 
physics beyond the SM. It provides a natural solution to the hierarchy problem 
and allows the unification of gauge couplings in the context of grand unified theory 
(GUT). Moreover, in the model with conserved R-parity, the lightest supersymmetric 
particle (LSP), typically the lightest neutralino, is stable and provides a dark matter 
candidate. 

In general SUSY models, the non-vanishing off-diagonal elements of the slepton 
mass violate the lepton-flavor conservation and then are strongly constrained by 
experiments ||. Since arbitrary sfermion masses generally lead to flavor-changing- 
neutral-current (FCNC) processes exceeding the experimental bounds, universality 
among generation of sfermion masses is usually imposed at some high energy scale. 
The most popular one is the minimal supergravity scenario (mSUGRA) ||, whose 
input parameters are the universal scalar mass m , the universal trilinear scalar 
coupling A , the gaugino mass mu 2 at the GUT scale and tan/3. 

Only small regions of mSUGRA parameter space have DM relic density within 



the observational bound fL0|| . This is because >m 1( / 2 is predicted for moderate 
rriQ ~ mi/2 and the LSP is almost pure Bino in this scenario. For moderate tan (3, only 
t-channel Ir exchange diagram gives important contribution to DM annihilation cross 
section (Tl[ . Therefore, the cross section is typically too small and the relic abundance 
exceeds the experimental upper bound. However, when relaxing the universality 
condition of the SUSY breaking parameters at the GUT scale, one can freely adjust 
the thermal relic density. 

One of the well-motivated relaxations is to allow soft scalar Higgs masses to be 



varied such as in the non-universal Higgs mass model (NUHM) JT2[ Jl3| [Tj] . 

Supersymmetric particles may be discovered in the early stages of LHC data 
collection. The squarks and gluino with masses below 1.5 TeV are expected to be 
found at the LHC for £ = 1 fb _1 at ^fs = 14 TeV [16||. By studying the kinematics of 
long cascade decays, one can extract information on masses of sparticles involved and 
more or less constrain the parameter space of the underlying theory. This is especially 
successful when sleptons are involved in the cascade decays. However, it was pointed 
out that, for a given set of experimental SUSY signatures at the LHC, more than 
one set of parameters could be present and ambiguities predominantly occur in the 



electroweak-ino sector. This is sometimes called the LHC Inverse Problem []17|. A 
similar result was obtained in [18|] when an exclusive likelihood map of SPSla point is 
studied by a weighted Markov chain technique. Even for the favorable SPSla point, 
we have multiple solutions. 

In this paper, we study the interplay between the LHC measurements of the pa- 
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rameters in electroweak (EW) sector and lepton-flavor violation. The LFV processes 
are important discovery channels of physics beyond the standard model. Right now, 
the main constraints come from the upper bounds of rare decay process searches 



BrQi -> ej) < 1.2 x KT 11 , (1.1) 
Br(r -> fij) < 4.5(6.8) x 1(T 8 , (1.2) 
Br(r -> ej) < 1.2(1.1) x 1(T 7 . (1.3) 

The bounds for tau LFV processes come from Belle (Babar). On-going MEG exper- 
iment |22| aims to push the sensitivity of fi — > e7 down by two orders of magnitude, 



and the SuperKEKB |23j and SuperB Factory [£4j] aim to improve the sensitivity of 
the r decay by one or two orders of magnitude in near future. 

If slepton is directly observed at the LHC, together with these advanced LFV 
measurements, we will be able to study the flavor sector of supersymmetry in detail. 



Several papers [15) |25| |26| recently point out that the ability of the LHC to measure 
the slepton mass differences will significantly constrain the possible scenarios. The 
further improvement may be obtained through the constraints on U(l) and SU(2) 
gaugino masses, M\ and M2, and fi parameter from the neutralino and chargino 
mass measurements at the LHC. Note that the inverse problem has an impact on 
the fj, — > e7 study. As will be discussed in detail later, the Br(p, — > ej) could 
be suppressed due to cancellation among different diagrams [p7| . The severity of 



the cancellation depends strongly on value of \i parameter and slepton masses. Ac- 
cordingly, to unravel the off-diagonal elements from the measurement of rare decay 
processes such as Br{p, — > ej), one needs to solve the inverse problem and determine 
the SUSY parameters precisely. 

In this paper, we take the non-universal Higgs mass model with positive \i with 
\i ~ m\/2 due to the DM relic abundance reason. In addition, we choose a mass hier- 
archy m^o < rrie R < m^o < rrie L < m^o so that both left- and right-handed sleptons 
appear in neutralino decay. Notably, we find mg R < m^o and m^ L < m^o should be 
satisfied when there is strong cancellation among the LFV diagrams involving the 
right-handed sleptons for the model points with correct DM relic abundance in the 
NUHM. For this parameters, many cascade decay modes involving sleptons can be 
observed at the LHC. Using the information, a solution of the EW parameters can 
be selected among the multiple solutions, and ambiguity in the MSSM parameters 
would be reduced significantly, leading better prediction to Br{ji — > ej) from the 
LFV slepton masses. 

The paper is organized as follows. In the next section, we introduce models we 
study. The sparticle mass determination is studied by explicit MC simulation in 
Section III for our sample model point. In Section IV, we discuss the ambiguity of 
the corresponding MSSM parameters and demonstrate that it can be removed by 
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studying decay branching ratios of Xi —>■ 21 + X, charge asymmetry of the decay 
distribution, and the rate into 2 hard jets + missing Et mode at the LHC The 
implication to LFV studies is discussed in Section V. We emphasize the importance 
of the slepton masses determination and compare the sensitivity to the off-diagonal 
elements of slepton mass matrices at the LHC and the other LFV searches such as 
MEG experiments. We also study the relation among slepton masses and \i parameter 
at the point of Br {fx — > e^y) cancellation. We show that the MSSM parameters would 
be studied very precisely at the LHC through the cascade decays involving sleptons. 
Finally, Section VI is devoted to discussions. 



2. Models 

The models under consideration in this paper consist of two important features. The 
first one is related to the scale of sparticle masses and their hierarchy. We chose the 
NUHM model as a representative. The second class involves issue of lepton-flavor- 
violation. Here, two models are selected: the MSSM with right-handed neutrinos 
and the MSSM with horizontal symmetry. Although models of interest have both 
NUHM and LFV, the detail descriptions of these features are given separately below. 

2.1 Non-universal Higgs Mass Model 

This model is motivated by the supersymmetric grand unified theories in which the 
Higgs fields do not belong to the same multiplet as the matter fields. Therefore, it 
is natural to expect that different multiplets would have different soft masses at the 
GUT scale. 

In mSUGRA scenario, for the correct pattern of EWSB, rn? Hu + /i 2 > — m|/2 
must be satisfied at weak scale. Since the SUSY breaking mass squared for the Higgs 
doublet Hu, m 2 Hu , is always driven to large negative value at weak scale by large top 
Yukawa coupling, is typically large. A smaller value can be easily obtained 
if one relaxes universality condition, especially mH v > rn Q at the GUT scale. Once 
~ mi/2, the LSP has larger Higgsino components. In that case, s-channel Z/h° 
exchange contributions and XiXx ~* WW contribution to the DM pair-annihilation 



cross section become more important [28| and the relic abundance would be small 



to be compatible with observational value. 

On the collider physics side, the reduction of would lead to a rich pattern 
of colored sparticle cascade decays at the LHC. Due to their considerable Wino 
components, Xa an d xt can be produced copiously in SU(2) doublet squark decay. 
Their successive cascade decays are the followings 



:2.d 



xt - 


-> vl - 


-» Xi , 


xl- 


*Ir- 


-xi, 


x2- 




- xi xl 
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which lead to clean, fruitful opposite sign same flavor (OSSF) dilepton events. The 
maximum values of dilepton invariant mass distributions depend on the cascade 
decay chains. If several endpoints can be identified and measured, it gives access to 
sparticle masses, or at least, their mass differences. 

2.2 LFV Models 

In the MSSM the low-energy LFV processes are induced by the off-diagonal terms 
in the slepton mass matrices, which depend on the origin of the SUSY breaking and 
physics beyond the MSSM. In this paper we discuss two models which predict LFV, 
i) MSSM with right-handed neutrinos and ii) MSSM with horizontal symmetries. 

Before going to those models, we first discuss /i — > process and an anomalous 
magnetic dipole moment a i} = (g — 2)i y The effective operators for Ij — > Z~7 and 
an are written as 



C eff = I^F^A^Pl + A*P R )k (2.2) 



where is a mass of the charged lepton lj. The coefficients Afj and A^j are functions 
of masses and mixings of sparticles inside the loop corrections. Their full formula 



can be found in [27]. The SUSY contributions to the anomalous magnetic dipole 



moment of the muon, a^ USY , is given by 



-ft •■"ft 

The branching ratio for the decay lj — > l~ r y is 

487r 3 a 



af SY = ml(A L 22 + Ag). (2.3) 



Br(lJ - 1,-7) = -^2-041 + \ A W) x Br(lJ - IJD^). (2.4) 



The experimental value of the anomalous magnetic moment of the muon mea- 
sured by the E821 experiment at Brookhaven with extremely high precision [2{| is 
given by 

a^ xp = 116592080(63) x 10 -11 . (2.5) 

According to the most recent calculations of the hadronic contribution based on the 
e + e~ data 1 301 , the difference between the SM prediction and the experimental value 



is 



Sa„ = a^ xp - af 1 = +302(88) x 10 -11 , (2.6) 

which corresponds to 3.4-a deviations. It was pointed out that this discrepancy 
tends to come from new physics contributions as it is unlikely to be explained by 
errors in the determination of the hadronic contributions pl| . Moreover, in the 



context of SUSY, the 5a^ anomaly, Eq. (|2.6|) , suggests light slepton-chargino sector. 
Accordingly, unless flavor mixings are very small, \i — > c-f may be detected at MEG 
experiment |22[ soon. 

We now explain two LFV models and discuss their phenomenology. 
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2.2.1 MSSM with Right-handed Neutrinos 

This model is motivated by the observations of non- vanishing neutrino masses and 
neutrino oscillation. In the model with three right-handed neutrinos, the superpo- 
tential for the lepton sector is given by 

W = iri-:;!,,!!,, + fjN^Hu + ±M?NtNf, (2.7) 

where L i; Ef and Nf represent chiral multiplets for left-handed lepton doublet, right- 
handed charged lepton singlet and right-handed neutrino singlet respectively, and Hjj 
and Hd for two Higgs doublets with opposite hypercharge. The neutrino mass matrix 
can be obtained by the seesaw mechanism and is given by 

m v = fjM-'Uihu) 2 (2.8) 

where (hu) is the vacuum expectation value of the neutral Higgs component h\j of 
the Hjj multiplet. 

In general, the Yukawa couplings fi and f u cannot be diagonalized simultane- 
ously. Therefore, even though one assumes that mass matrices for left-handed and 
right-handed sleptons, m~ and m|, are proportional to the unit matrix at the GUT 
scale, the LFV masses of m| will always be generated via the renormalization group 
evolution, 

+2(/>4/„ + rhljtU + 4A0{] , (2.9) 

by the non-vanishing off-diagonal elements of f v and of the trilinear soft parameter 
A v WA - The large mixing angles observed in the atmospheric and solar neutrino 
oscillation experiments enhance the LFV masses of m~, and the experimental bounds 
on Br(fi — > cy) and Brir — > /ry) give stringent constraints on the model f33fl . 

It has been studied that if the LFV is generated only in the left-handed slepton 
sector, Br(fj J — > ej) will be strongly correlated with (a^ USY ) 2 [0. This is because the 
chargino-sneutrino diagram dominates over other contributions in both observables. 
By taking a common mass for all sparticles, the current experimental bounds on 
a^ USY (= <5a^) and Br (ft — > ej) put a stringent constraint on the left-handed LFV 
slepton mass term, 

<2xl0-'f '< SY| ,V7 Br( ^^y. (2.10) 
^3.02 x 10" 9 J \1.2x 10- 11 J V ; 

On the other hand, we still have a room for sizable mass difference between 
the left-handed selectron and other sleptons, which may be measured at the LHC. 



2 

m SUSY 
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Without LFV, basically the "left-handed" tau slepton mass is deviated from other 
left-handed slepton masses due to its Yukawa coupling effect to the renormalization 
group evolution or the left-right mixing term, while mass difference mg L — m^ L is 
negligibly small. However, in the case of nonzero LFV, mass splittings could be 
larger and would be detectable at the LHC. 

2.2.2 MSSM with Horizontal Symmetry 

Horizontal symmetries are introduced to derive the hierarchical structure in the 



Yukawa coupling constants |35|] and also to suppress the off-diagonal terms in the 



sfermion mass matrices [315] [17]. 

In this paper, we studied the MSSM with horizontal symmetry, which is origi- 
nated from one of models in []38[] . The original model was a supersymmetric model 



in which the SUSY breaking is mediated to the MSSM by gravity interactions. The 
gravity-mediated contribution is assumed to be subject to an approximate horizontal 
symmetry. More specifically, consider a horizontal U(l) x U(l) symmetry where each 
£7(1) is explicitly broken by a scalar singlet spurion carrying the corresponding charge 
— 1. The sizes of both U(l) breakings are assumed to be equal and then parametrized 
by a single parameter e ~ |V^ S | ~ 0.2. The horizontal charge assignment for lepton 
sector is as follows: 

^(4,0), L 2 (2,2), L 3 (0,4); 
#l(1,0), £ 2 (l,-2), £ 3 (0, -3). (2.11) 

Using above horizontal symmetry, we parameterize the left- and right-handed 
slepton mass matrices at the GUT scale as 



m- 



2 -ml + xmlX'j, M? = m 2 + xm 2 X r s , (2.12) 



where x is the ratio between the flavor-independent and dependent contributions. 
The structure of matrices X' L and X' R can be determined uniquely by selection rules 
by the horizontal symmetry and have the following forms: 



(2.13) 



Here, we neglect the flavor-dependent contribution to the flavor-diagonal mass terms, 





though they are also allowed by the symmetry with size ~ xm\. It is argued in J38 



that even x ^ 1 is allowed from phenomenological constraints under the horizontal 
symmetry when m < 380 GeV, my 2 < 160 GeV, and 5 < tan/3 < 15. 

To discuss the fi —>■ ej phenomenology of this model, first of all, it should be 
remembered that the effective /i — > operator, Eq. ( |2.2|) , must flip the chirality 
and change the flavor of the external leptons. Since the left-handed LFV masses 
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are suppressed in this model, the dominant LFV contribution comes from the right- 
handed sector. In the mass insertion approximation, there are two dominant one-loop 
diagrams when M\\x tan j3 is large |[39| . The former one has chirality flip on an internal 
line via left-right slepton mixing and its contribution is 

al i lay 1 H)^M 1 (A M + //tan/3) 
12l5 °" 2 47rm? m 2 mf ' 1 j 

In the latter diagram, the chirality is flipped at one vertex via Yukawa coupling. 
The contribution is proportional to the Higgsino components, which come from the 
Bino-Higgsino mixing of neutralino. The amplitude is given by 

al I a Y 1 ( m f fl )e Mi/i tan/3 

A2l(BO-5o mixing) - 47r _ 2 _ 2 n h{- 2 )- {^) 

l R IR IR lR 

In the above equations, m~ (m~ ) stands for an averaged right-handed (left-handed) 
slepton mass, and (m~ is a (/t, e) component of the right-handed charged slepton 
mass squared matrix. The kinematic function f 1 is given by 

8-llx + 4x 2 -x 3 + 2(2 + x)logx , n . 

= 2(1^)4 ( 2 - 16 ) 

which is a positive-definite, decreasing function of x. 

A key point is that the relative sign of the two amplitudes is negative and then 
cancellation between diagrams can occur significantly. For = 0, a severe cancel- 
lation happens when 

1 1 u 2 
- — o o-A(-V) ~ 0. (2.17) 

2m? m? v m 2 ; v ; 

li lR IR 

This cancellation occurs when m 2 ~ u 2 in the limit of u 2 /m 2 ^> 1. Since dominant 

i>l ' Ir 

contributions for fi —> ej and a^ USY are now different, then their correlation becomes 
much weaker than in the case that the left-handed sleptons have LFV masses. 

Next, we consider the processes /t — > eee and \x — > e conversion in the nuclei. 
Their formulae are also given in P? |. See also Ref. for precise evaluation of 



/j — * e conversion in the nuclei. These two observables also suffer from the partial 
cancellation [39 1 . The /a — > eee receives dominant contribution from the penguin- 



type diagrams which are enhanced at large tan/3 region in a similar manner to 
/i — > e7 so that their behaviors are alike; the subdominant box-type contribution 
just helps lifting up the depth of the cancellation valley. On the other hand, the 
jj, — > e conversion rate behaves rather differently and plays a complementary role in 
resolving the cancellation point. Thus, the correlation among the LFV processes are 
also weaker than the case that the left-handed sleptons have LFV masses. 
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m Q 


100 


mi/2 


300 




100 




380 


tan p 


10 




271.33 


Mi 


122.49 


M 2 


230.89 


m- g 


719.67 
















665.19 




648.85 


m d T . 


670.29 


m 4 


642.47 


m w 


600.91 


m h 


638.72 




462.35 


m t2 


655.20 


m e L 


239.62 


m e R 


130.38 




128.07 




238.89 


VTLy 


224.37 


mc T 


222.16 


Xi 


196.30 




321.62 




114.70 




197.82 


Xa 


278.87 


A4 


323.23 




111.22 




350.05 


m A 


347.31 


m H + 


358.56 



Table 1: Relevant parameters and sparticle masses in GeV for the parameter point A. 



3. Monte Carlo Study 

In this section, we study the leptonic SUSY signals at the LHC in the non-universal 
Higgs mass scenario. For our analysis, we use ISAJET v7.75 |^1[ to calculate the 
sparticle spectrum and IsaReD P2"| , which is part of the IsaTools package, to evaluate 
the dark matter relic density. We generated 5 x 10 6 events by HERWIG 6.5 ^3[]; this 
corresponds to about 300 fb -1 of integrated luminosity. The AcerDet package [|4[] is 
employed to simulate the detector response. 



3.1 Model Point 

As mentioned earlier, we chose m>H D = mo ^ mn u at the GUT scale. We are inter- 
ested in the case that both left- and right-handed sleptons can be directly produced 
via neutralino and chargino cascade decays. Then slepton masses should be light and 
a relatively small value of m had been chosen. Furthermore, we chose a moderate 
value for m\/2- The relevant parameters and sparticle masses for our studied point 
A are listed in Table § We took m = 100 GeV,m 1/2 = 300GeV,m^ = 380 GeV. 
This leads to /i = 271 GeV, mg L = 240 GeV, and m~ eR = 130 GeV. We will see 
that this corresponds to the point where strong cancellation among contributions to 
Br(fi — > ej) occurs. 

For this choice of parameters, the value of fi substantially reduces from = 
397.30 GeV for mSUGRA case with inputs: m = 100 GeV, = 300 GeV, and 
tan/? = 10. This increases the Higgsino components of the LSP as the mixing 
matrix elements (N s , N w , N Sd ,N Bv ) = (-0.96, 0.09, 0.22, -0.11). The dark matter 
relic density is evaluated to be 0.1179 which is consistent with the combined results 
from WMAP and SDSS @§, tt DM h 2 = 0.111+°;°^ (2a). 

The relevant sparticle decay branching ratios are compared with mSUGRA case 
in Table [| In mSUGRA scenario, Xi 1S almost pure Bino, and x\ an d xf are 
Wino-like. Then q L decays substantially into x\ an d X\- As the decay into is 
kinematically forbidden for this parameter choice, they essentially decay into left- 
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A 


mSUGRA 




A 


mSUGRA 


ul -> xl 


25.4 


31.6 


xl -> e L 


2.3 


1.0 


ul -> Xs 


0.2 


0.1 


xl £r 


0.7 


0.3 


-> X4 


7.8 


1.3 


xl -> e fi 


13.5 


3.1 


«l -> Xi 


53.2 


64.4 


xt 


7.4 


2.1 


ml -> xt 


13.1 


1.8 




88.8 


63.3 


d L -> X° 


21.8 


30.7 


Xl" 


26.1 






0.4 


0.1 


ei -> xl 


26.5 


5.9 


-> X4 


9.7 


1.8 


zl -> Xi 


35.9 


84.0 


rfi -» xr 


42.1 


60.4 


e/? -> Xi 


100 


100 




23.1 


4.8 









Table 2: Relevant branching ratios in % for our sample point A compared with mSUGRA point. 



handed component of f\ with large branching and small Bino component of xl could 
decay into e# producing the well-known edge of dilepton invariant mass distribution 
On the other hand, heavier inos are Higgsino-like and would not be produced 
at the LHC This will then forbid cl to show up at the LHC also. 

For point A, when \x is smaller, heavier inos have larger Wino components so that 
their productions in the left-handed squark decay become significant. Table [| shows 
that Br(qi — > xl/xt)i along with Br(xl — * &l) and Br(xt ~* ^l), are enhanced 
considerably. On the contrary, the enhancement in Br(\l ~^ &r) well demonstrates 
the increase of Bino component in xl- Summarily, with small a ~ mi/2, a number 
of sparticles could show up through various decay patterns, Eq. (|2"Tl"|), at the LHC 
experiment. Note that X3 is almost pure Higgsino and then is not produced neither 
at point A nor mSUGRA point. 

3.2 Two-lepton channel 

In analyzing SUSY signals at the LHC, we put an emphasis on lepton channels 
when lepton presumably means electron or muon. In this section, we focus on the 
celebrated dilepton invariant mass distribution and model independent constraints 
on the sparticle masses. The other signatures will be discussed in the next section. 
The dominant neutralino and chargino cascade decay processes which lead to OSSF 
lepton pair in the final states and the corresponding expected kinematics endpoints 
are listed in Table |3] 



The events are selected by the following criteria 

• an OSSF dilepton pair where both leptons have p l T > 10 GeV and \r)\ < 2.5, 

• more than 4 jets with p 3 T1 > 100 GeV,p^ 2 3 4 > 50 GeV, 

• M cS = pr A + p T ,2 + Pt,s + Vta + tf T > 400 GeV, 
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decay mode 


m m 


^min 
m jU 


m jl 


m min 


m u 


(1) Ql -> xl -* -> Xi 


517.4 


209.1 


477.6 


272.7 


70.7 


(2) Ql * X.2 * * Xi 


461.2 


202.8 


417.2 


253.8 


111.6 


(3) <?l -> X°4 -> e L -> X° 


459.8 


223.4 


390.2 


285.7 


122.4 


(4) <7l -> X° e fl -> X? 


550.3 


217.1 


532.0 


249.6 


140.6 


(5) q L -> X° g L -> X? 


550.5 


296.5 


510.4 


383.6 


190.5 



Table 3: Endpoints of invariant mass distributions in GeV for various decay processes. 




Figure 1: Dilepton invariant mass distribution (in GeV). 



• $ T > max(100,0.2M eff ). 

The dilepton invariant mass distributions are shown in Figure [I]. In these plots, 
the subtraction of opposite sign opposite flavor (OSOF) dilepton distribution is uti- 
lized to reduce SUSY backgrounds. The endpoint position from c[l — > Xi ~^ I ~* Xj 
is given by analytical formula p5f 



(m% — m~)(m~ — m% ) 
m u -y -2 (3-1) 

From the dilepton invariant mass distribution, edges of all decay modes except 
(3) in Table |3| are visible. The distribution of decay mode (3) is small because 
basically x° decays further through the golden mode (1) and results in four-lepton 
final states discussed below. 

3.3 Four-lepton channel 

Thus, we turn to consider four-lepton events. If these leptons really come from the 
cascade decay 

Xl - e L - *2 - e R - xi (3.2) 

they must form two OSSF pairs, and an invariant mass of one pair must be below 
the endpoint of mode (1) while another be lower than that of mode (3). This mode 
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Distribution 


Expected endpoint 


Fitted endpoint 


m U (i) 


70.78 


70.67 ± 0.01 


m U (2) 


111.61 


111.70 ± 0.21 


mi(3) 


122.41 


121.93 ± 0.43 


m U (4) 


140.62 


140.64 ± 0.26 


m U (5) 


190.46 


191.05 ± 0.28 



Table 4: Fitted endpoints in GeV 



is very useful in the two senses. Basically, ~ ^x2' m xt ~ 



m^o and ~ m gi so 



that the endpoints from modes (2) and (3) are close. If endpoints of other modes are 
much different from these two, by looking at from four-lepton events, one can 

easily pinpoint which edge from dilepton events is the edge of mode (2). Moreover, 
we can use as a cross check for m^o and mg L obtained from edges of mode (4) 

and (5). 

We now consider the possibility to identify the decay mode (3) whose daughter 
X2 subsequently decays via mode (1) producing four leptons in the final states. The 
following cuts are applied to select events: 

• exactly two pairs of OSSF leptons where each lepton has p l T > 10 GeV and 

M < 2.5, 

• four leptons must be composed of exactly one e + e~ pair and one pair. 



Entries 

Mam 

RMS 



IS15 
109.0 
17.55 




120 140 160 

win (41 events) 



Figure 2: Dilepton invariant mass 
mu distribution (in GeV) from four- 
lepton events with the requirement 
on the other lepton pair mini < 
70.7 GeV. 



The m ee and are then calculated for 
each events. After requiring that an invariant 
mass of one OSSF lepton pair must less than 
70.7 GeV, distribution of the other pair is shown 
in Figure §. Now the edge around 122 GeV shows 
up confirming that events are truly from the de- 
cay in Eq. (|3.2| ). For a consistency check, we 
also plotted the trilepton and four-lepton invari- 
ant mass distributions from events which pass 
the above cuts; the plots are shown in Appendix 
A. The fitted values of various endpoints are ob- 
tained by using linear fit function 



f(M) 



AM + B, 
0, 



< M < M™ ax 
M > M/f ax 



(3.3) 



smeared with a Gaussian and are listed in Table 
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m- qL 


Solution 1 


114.70 


130.38 


197.82 


665.19 


Solution 2 


178.09 


244.03 


265.06 


747.22 



Table 5: Sparticle mass ambiguities in GeV 



3.4 Mass Parameter Determination 

Once kinematical endpoints have been measured experimentally, relevant sparticle 
masses may be extracted. It is expected that m^, m^o, rn~ lR) and m^o can be recon- 
structed from m™ ax , m™ x , mj 1 ^ and m™ n j|7|. Furthermore, if higher endpoints 
are visible, one would then be able to identify the decay of heavier neutralino and 
then resolve masses of xl an d II- 

For a given set of m^ ax , m^ ax , m , Jl ax and m™ n measurements, in principle, there 



may be different corresponding sets of sparticle masses [ 4"8"|| . For the set of endpoints 
(1) in Table |3], there are two possible solutions as listed in Table |^. The second 
solution can be discarded by measuring rn^f 1 , which differs by around 20 GeV. 
Alternatively, solution 2 does not have solution for rrig L which satisfies rn™?fi and 
m™ 5 x simultaneously. 

In order to estimate error for each SUSY particle mass, we generated a set of 
random numbers corresponding to a set of masses {m^o, m^o, m^o, m^, mg L }. We 
then calculated a set of measurable quantities {m^ n , m^ x , mj[ m , m™ x , m™^ 4 
and their chi-squared which is defined as 

. 9 (Nominal value — Measured value) 2 .„ JS 

Ax 2 = V ^ 2 3.4 

erf 

all observables ' 

In the above definition, <7j includes both systematical error and statistical error. 
We employed systematical errors as listed in Table |6| ||49|| . In our analysis, however, we 
used only statistical errors for rn^^ 4 5 % since they dominate over systematical errors 
(see Table |]). Moreover, we fitted endpoints for m^ m , m^ ax , m™ n , m^ ax distribu- 
tions and found that even if they have good statistics, their fitted values curiously 
differ from central values, especially m™^ x . We then just used their systematical 
errors in calculating A% 2 . 

In addition, the negligibly smallness of both systematical and statistical errors for 
implies that will be measured rather precisely. Therefore, in the analysis, 
we took m™/"? as an input and obtained rrie R for given values of {m^o, m^o}. Results 
of l-o" error estimation for sparticle mass differences are shown in Table [7[ 



4. Flipping Solutions 
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Sparticle Mass 


Central value 


Estimated error 


A.1 


114.70 


+6.7 
-6.3 


rrie R - m^o 


15.68 


+0.45 
-0.49 


m^.o — m^,o 

A.9 A] 


83.12 


+0.75 
-0.62 


mg L — m^o 


124.92 


+0.65 
-0.65 


m^.o — m^o 


208.53 


+0.77 
-0.64 




551.19 


+4.64 
-4.47 



Table 7: Central values and 1-a error estimation of relevant sparticle masses in GeV. 



Distribution 


Systematical error (GeV) 


m max 


0.08 


m max 


4.3 


m ill 


2.0 


m max 


3.8 


m mm 


3.0 



Table 6: Estimated systematical errors for different 
endpoints (in GeV). 



Besides the one shown in Table 
H, there is another kind of ambigu- 
ities which is related to parameter 
point identification especially when 
the ordering of fi, Mi and M 2 is shuf- 
fled. We now introduce flipping so- 
lutions. Flipping solutions are so- 
lutions among which masses of rel- 
evant sparticles (masses of left-handed 
squark, left- and right-handed slep- 
tons, and three neutralinos which have significant gaugino component) are the same 
but ordering of Mi, M 2 and \x parameters are different. The mass degeneracy among 
solutions results in the same endpoint positions from those sparticle decay. There- 
fore, only endpoint measurement is not enough to distinguish these solutions. By 
reminding that M\ < M 2 < /i for point A, we illustrate another two flipping so- 
lutions: Mi < < M 2 (point A2) and fi < Mi < M 2 (point A3); their relevant 
parameters are listed in Table ||. 

We fixed three masses of neutralinos which have significant gaugino component 
equal, as they are frequently produced from squark decays and would be measured 
rather precisely. For points A and A2, such three neutralino states are x? ; X 2 and xSL 
so almost all masses except m^o for both points are degenerate (m^o| P oi n tA = 279 GeV 



and m^o | point A 2 = 234 GeV). For point A3, however, they are x?>X3 an d xl instead. 
Therefore, we fixed m^o = m^o|p 0int A and m^o = 155 GeV. On the other hand, X3 
for points A and A2 and x 2 for point A3 are nearly pure Higgsino, and they are not 
produced from squark decays, so that the experimental constraint would be weaker. 

In the followings, we discuss and compare some phenomenological signatures of 
all flipping solutions. We will show that different properties of neutralinos essentially 
lead to discrepancy in both collider and low-energy LFV observables. 

4.1 Branching Ratios 

First of all, we show the OSSF dilepton distributions in Figure 0. The total number 
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A 


A2 


A3 




A 


A2 


A3 




271.33 


226.06 


146.21 










M 1 


122.49 


125.66 


187.46 


M 2 


230.89 


272.86 


291.47 


rriu L 


665.19 


665.19 


665.19 


man 


648.85 


649.84 


657.09 


m dj, 


670.29 


670.29 


670.29 




642.47 


641.67 


643.68 


m iL 


239.62 


239.62 


239.62 


m e R 


130.38 


130.38 


130.38 


m f2 


238.89 


270.32 


292.60 




128.07 


123.82 


173.80 




224.37 


224.03 


223.92 


m~g 


719.67 


719.35 


720.81 


m xt 


196.30 


193.34 


133.27 




321.62 


320.88 


319.79 


Xi 


114.70 


114.70 


114.70 


X?, 


197.82 


197.82 


155.17 




278.87 


234.28 


197.82 


m x° 4 


323.23 


323.23 


323.23 



Table 8: Relevant parameters and sparticle masses in GeV for points A, A2 and A3 respectively. 





A 


A2 




A 


A2 




A 


A2 


u L -> xl 


25.4 


15.4 


u R -> X° 


4.1 


9.6 


xl -> e R 


13.5 


14.1 


ml -» xl 


7.8 


17.8 


-> X° 


4.1 


9.6 


xt ^ 


7.4 


12.9 


ul -> xt 


13.1 


33.2 


x2-^ + x° 


0.2 


0.2 


x+ - Z + *+ 


16.1 


13.4 


dL -> X° 


21.8 


10.4 




1.4 


0.3 


^ -> xi" 


26.1 


20.4 


g?l -> Xa 


9.7 


20.5 


X° g L 


2.3 


4.2 


ei, -> X° 


26.5 


25.6 


dL X 2 


23.1 


45.2 


X4 -> g R 


0.7 


0.3 


ei, -> Xi 


35.9 


47.8 



Table 9: Comparison between relevant branching ratios in % for points A and A2. 



of generated events is 5 x 10 6 for all points. This corresponds to about 300 fb _1 
for points A and A2 and 250 fb _1 for point A3. The difference between distribution 
for points A and A2 is conspicuous. Since xt an d Xa are Wino-like at point A2, 
Br(q L — > xt/xl) an d Br(xt/xl Il/^l) are enhanced (see Table |9|). These result 
in the Z peak and the height of edges around 110 GeV (from xt ~ > ~^ Xi) an d 
190 GeV (from x\ ~> &l —> Xi) f° r A2 at least twice the size of that for point A. 

The enhancement of the branching ratios Br(q~L — > xt) an d Br(xt — > v£) can 
be understood analytically. Consider the x^-lepton-slepton interaction 

C oc g$m(j) R (v*x 2 P L e + ePgjfev) + g sin <f) L (e* L x 2 P L v + vP R xt^L)- (4-1) 

Here Pl,r stands for left-, right-handed projection operator and the mixing angles 
4>l and 4>r_ are given (with two- fold ambiguity) by 

2v / 2m W /(M 2 cos f3 + fj, sin/?) 
tan 2(t)L = Ml-^-2mi w cos2p ' ^ 

tan20 fl = 2 ^ (M2Sin/j 2 +/iCOS/3) . (4.3) 
VR Mf- Lt 2 + 2m 2 w cos 2(3 V ; 
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The interaction x.2 ~ v — e and, analogously, X2 — u — d are proportional to sin 0# 
which becomes larger when M 2 and /i are inverted so that \i < M 2 , for moderate 
tan (3. 

Point A3 can also be distinguished from the others. The important key is that, for 

point A3, mass difference m-+ —m~+ is close to m^o —m^o, in addition to m^o —m^o. 

Hence the dislocation of rn^ ax (xt ~* &l Xi) w i n serve as an indicator of flipping 

solution. In Figure |3], the highest endpoint around 180 GeV for point A3 actually 

comes from the decay xt ~ ¥ &l — > Xii while the distribution for xl ~ ¥ &l — ► Xi 
(jYfjnax ^ -^gg Q e y) j s muc h smaller. 

In addition, the study of four-lepton events in analogous to subsection [T3] re- 
veals the endpoint from the decay chain xt ~ * &L —> X3 {rriu' ax ~ 120 GeV) where 
X3 subsequently undergoes the cascade decay X3 ~ * &R ~ * Xi- By utilizing all end- 
point information from four-lepton events, X4 an d &l masses can be extracted and 
m™ x {x\ —> &l — > Xi) ~ 190 GeV can be calculated. To this end, one can confirm 
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205.64 


M 1 


253.75 


M 2 


414.91 








747.54 




696.83 


m h 


752.02 




680.72 


m iL 


389.37 


m e R 


243.92 




387.60 


m fl 


233.51 


mo 


380.13 


m~ g 


766.77 


m-+ 
Xi 


198.64 




435.96 




179.04 


X.9. 


212.99 


A3 


265.21 


m x° 4 


440.05 



Table 10: Relevant parameters and sparticle masses in GeV for point A4. 



that the win endpoint at 180 GeV does not come from x° decay. 

A very striking feature of its double peak distribution below the Z peak in Figure 
|3] seems to make fi < Mi < M 2 case very easily distinguishable from other cases. 
However, this actually depends on our model parameter choice, i.e. m 



'<--R 



e R - 



The distribution receives contribution to the outer edge (around 70 GeV) 
from the Bino-like x® cascade decay x® ~~ * Ir ~ > Xi an( i to the inner edge (around 
40 GeV) from the Higgsino-Bino mixed x° cascade decay Xz ~* Ir ~* Xv Consider 
the case that our true model point is the solution 2 in Table [5[ Its fi < Mi < M 2 
flipping solution (point A4) has m^o < m^ R and an unusual double peak distribution 
disappears. We show mass spectrum of the point in Table |l^ and ma distribution in 
Figure £|. 

4.2 2 jets + Ij/rp Signature 

The distinctive Higgsino-like property of Xi f° r 
very small \i solution makes it possible to dis- 
tinguish one from intermediate-to-large \x solu- 
tions. More specifically, it is clearly presented in 
the differences in the decay branching ratio 



Entries S67957 
Mean 49. 74 

RMS 21.49 



Br(qi 
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Figure 4: OSSF dilepton invariant 

mass distribution for point A4. 

Pt jets + large lf T and no isolated 

number of jets. We select events by 

(cl) more than 2 jets with \r]i\ < 

200 GeV, 



point A 
point A2 
point A3 

When both pair-produced right-handed squarks 
decay directly into Xi> their signature is two high 
lepton in the final states. Here, we first count 
employing the following cuts: 

2.5 and transverse momentum p-r^ > 250 GeV,pT,2 > 



(c2) no isolated lepton with px greater than 10 GeV, 
(c3) no tagged b-jet, 
(c4) tf T > 250 GeV. 
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im> 

1600 
1400 
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800 
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Figure 5: The number distributions of jet with px > 50GeV after the event cuts for point 
A and A3, respectively. The dashed line represents the distribution for pp — » quqn — ► XiXi- 
Notice that the vertical axis of the left-handed figure is multiplied by 100. 




mn (npo =2) m-n (n j50 =2) 

Figure 6: tot2 distribution for nj$o = 2. The left and right panels are for point A and A3 
respectively. The dashed line represents the distribution for pp — ► q~nq~R — > XiXi- 

The number distribution of jet with pt > 50 GeV, n^o, after cuts is shown in 
Figure || Note that the distribution for point A2 is very similar to that for point A 
and therefore not shown here. There are obviously two differences between these two 
plots. Firstly, high jet-px cut (cl) and lepton cut (c2) reduce the number of events for 
point A3 so substantially that it is an order of magnitude smaller. Secondly, point A3 
has smaller fraction of events for n^o = 2 and larger fraction when rij^o > 4. These 
all indicate that squark tends to decay via a longer cascade decay which produces 
leptons or several softer jets in the final states. It further implies that squarks would 
have considerable decay branching ratio into heavier inos and the LSP has small 
gaugino components, i.e. Higgsino-like. To this end, one can infer that point A3 has 
small [i < Mi, M 2 . 

In the same plots, the dashed line represents distribution of pp — > qnqu — > XiXi 
events. Even though there must have only two jets at parton level, the initial state 
radiation can contribute to the high px third and forth jets. The difference between 
the solid and dashed line receives contributions mostly from left-handed squark or 
gluino (associated) production. 

The method of counting jet numbers, however, is not so reliable. We need more 



Entries 5050X4 
Mean 3.650 18000 
RMS 1.112 

16000 
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X?) 



concrete evidence to confirm the conclusions about small Br(q 
small value of \i. This can be achieved by the so-called mri method |52 
variable is defined in event-by-event basis as 



and hence 
. The rriT2 



raraCPr ' Pt> P^T 1SS i m test) = a min max |m T (p^, j^; m test ), m T (p^ 2 , pff ; m tcst Xi4) 

where m test is a test mass and the minimization is performed over all possible splitting 
$t + P^t = Pr"^- The transverse mass my is defined as 

™t(Pt> Pt; m test) = + m 2 test + 2 - p 3 T ■ . (4.5) 

The rriT2 variable has the property 

m T 2(m test = m^o) < rrig R , (4.6) 

for pp — > g_Rgi? — > events, so that one can determine from the endpoint 
measurement of rriT2 distribution. 

Figure |6] shows m T2 distributions for rij 50 = 2. Here, m test is taken to be the 
nominal value of m^o. Clearly, both distributions with n^o = 2 have endpoint 
about the correct qu masses. However, the distribution for point A has sharp edge 
and events near the edge come mostly from the true cascade while that for A3 is 
contaminated mostly by contribution from qi production. 

In order to estimate a statistical significance of signal over SM background, we 
adopted a set of event selection cuts from [ |53[| : 



• ^ r > max(200 GeV, 0.25M e// ) and M eff > 500 GeV, 

• two jets with p T > max(200 GeV, 0.25M e //), \r]\ < 1 and AR > 1, 

• no additional jet with p T > min(200 GeV, 0.15M e ff), 

• no isolated leptons and no tagged b-jets, 

• transverse sphericity St > 0.2. 

It should be noted here that the above two hardest pt jet cut, pr > max(200 GeV, 0.25M e //), 
distorts shape of m T2 distribution. This is because the event near m T2 endpoint cor- 
responds to the configuration where two jets go in the same direction with ptj ~ 
0.25M e //. Therefore, this cut kills significant numbers of events near the endpoint. 

In Figure [j], we just show how rriT2 distribution is distorted. On the left panel, 
rrtT2 distribution for point A is shown when the above cuts are applied at ~ 30 fb^ 1 
of integrated luminosity. It is an M e ff dependence of the cut that selects events with 
rather high pt jets. If this cut is relaxed to be px,i > 250 GeV,pr,2 > 200 GeV, 
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mn inn 

Figure 7: rriT2 distribution of point A at 30 fb^ 1 of integrated luminosity. In the left 
panel, selection cuts are taken from []53[| . The two hardest px jet cut is different for the 
right panel as described in the text. Again, the dashed line represents the distribution for 
PP -> WW -> XiXv 



Point 


No. of Signal 


No. of SM Background 


S/B 5M 


S/ V^SM 


A 


1341 


180 


7.5 


100.0 


A3 


133 


180 


0.7 


9.9 



Table 11: The number of signal and SM background events passing the selection cuts described 
in the text, signal-to-SM background ratio, and signal statistical significance for points A and A3 
at 10 fb~ x of integrated luminosity. 



in the right panel, the peak of rriT2 distribution leans more toward its edge as it is 
supposed to be. Moreover, the number of events passing the cuts are about three 
times larger with increasing fraction of pp — > q~nq~R — > XiXi events near endpoint. 

The signal-to-background ratio and signal statistical significance of events pass- 
ing above cuts for 10 fb~ l of integrated luminosity are given in Table |TT| for points A 
and A3. The background statistics is taken from f53|| , page 1630. The signal statistic 
of point A3 is so poor that its distribution is submerged in SM background. 
This makes it even easier to distinguish between two points. 

Summarily, by investigating the rriT2 distribution and branching ratio of heavier 
ino decay, the region to which the \i parameter belongs can be ascertained. 



4.3 Charge Asymmetry 

In this subsection, we will illustrate that a charge asymmetry between the jl + and 
jl~ invariant mass distributions can help lift degeneracy between points A and A2. 
This method was firstly proposed in |)IJ where the decay q~i — > X2 ~ ¥ h * s studied in 



the case that slepton is purely right-handed; later the left-right mixing was taken into 
account in (J5TJ. In our study, we assume that left-right slepton mixing is negligible. 
We now briefly explain the method following prescription in |50 . 

Consider the cascade decay q a — > g%° — > ql nea 4/3 where l neaT denotes lepton from 
X° decay and a, (3 = L, R. The ql nea r invariant mass is given by a simple kinematics 
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expression 



(m q i near ) 2 = (mZ:J 2 s 2 e/2 (4.7) 

where sq/2 = sin(0/2) and 9 is the angle between quark and lepton momenta in xl 
rest frame. Due to the chirality structure of quark-squark-neutralino coupling, Xi is 
polarized. Its polarization alters the angular distribution of its daughter lepton and 
hence the angular distribution of ql n ear invariant mass. Mathematically, the prob- 
ability density function receives extra angular-dependent factors from the chirality 
projector. For the case a ^ j3, the probability density is given by 



dP = \ 4s^ /2 for ql+ ear or ql near 

dse/2 \ 4s0 /2 (l - s 2 e/2 ) for ql~ ear or $1+^ 

If a — f3, the density function for ql~ ear and qlnear are inverted. Based on the fact 
that valence quarks have harder PDFs than sea quarks, squarks would be produced 
via qg production more than anti-squarks. In the following, we therefore assume that 
high-pr jet represents quark rather than anti-quark. 
Define charge asymmetry 



A = — where s = <icr(m„z±). (4.9) 

s + + s~ 



in'" 1 ; 



ql 



We obtain 



For the celebrated decay mode qi — > X2 ~ * Iri ^near with positive charge rather 
than negative charge favors to go to the opposite direction to the quark jet and 
constitutes events near endpoint. The charge asymmetry for this case is then positive. 

In our study, once we know that the LSP is Bino-like by mn method and that 
qL cascades via long decay chains, this method can be used to probe the property 
of X2- Figure ^| shows N(jhl + ) — N(jhl~) as a function of rrij h i for points A and A2 
at ~ 60 fb^ 1 of integrated luminosity. In these plots, we use only jh, one of the two 
highest pt jets which gives higher rriju value. Moreover, rrij h i is calculated only for 
lepton with higher p T because m^o — m e - fl > m e - fl — m^o and lepton from X2 decays 
tends to have larger px- In addition, we require mu < 70.7 GeV and rriju < m™^* . 
For point A (left), the charge asymmetry is clearly positive which indicates the 
q~L —>■ X2 ~ y decay. On the right panel, the asymmetry is also positive but less 
prominent. At rrij h i > 300 GeV, the distribution is relatively flat compared to the 
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Figure 8: The difference between m,j k i+ and rrij h i- distributions as a function of rrij h i for 
points A and A2. 

left panel. This can be interpreted as the contribution from q R decay. The actual 
ratio of Br(q R —> x^/Br^qi — > X2) * s 

Br(q R -> x§) _ f 0.16 point A 

5r(g L -» x°) ~ I 0.67 point A2 

The ratios differ by factor of four: a factor two from the reduction of Br(q~L —>■ X2) 
and another from the enhancement of Br(q R X2) ( see Table 0). This makes the 
contribution from q R for point A2 is comparable to q^s contribution and really cause 
the decrease in the charge asymmetry. Therefore the flatness of charge asymmetry 
implies that X2 have smaller Wino component and it is persuasive to conclude that 
Mi < fi < M 2 for this point. 

In Figure we show the dependence of R = Br{u R — > X2)/Br(u L — »■ %°) 
on ///Mi and M2/M1. Here, we took parameter point A and allow /x and M2 to 
deviate within 20 percent. The ratio R shows little dependence on MijM\ while 
it ranges from below 0.1 to 0.6 when fi/Mi ratio varies from 1.5 to 2.7. Due to 
the form of neutralino mass matrix, Wino mixes with Higgsinos easier than with 
Bino. Therefore, as a probe into the reduction of Wino component of the ratio R 
depends on fi/Mi stronger than on M2/M1. When fi/Mi is smaller, mixing between 
neutralino states is larger and the charge asymmetry receives more q R contribution 
so that the distribution becomes flat. 

Now we give comments on other possible nipping solutions. The first possibility 
is a solution when Mi < M2,// and mg R and rrte L are flipped. This case can be 
easily distinguished by utilizing a charge asymmetry to confirm the chirality of a 
daughter slepton in the decay qi — > x® ~~ h- However, if not only mg R and mg L , 
but also rriq R and m^, and Mi and M 2 are inverted with M 2 < M 1; /i, events near 
m™ ax ~ 71 GeV comes from q.n II Xi instead and charge asymmetry of 

the original and this solutions must be the same. Nevertheless, whenever Mi and 
M2 are flipped, the chargino sector will be affected. For example, for M2 < Mi < fi, 
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M 2 /M 1 

2.14 2.16 2.18 2.2 2.22 2.24 2.26 2.28 
0.6 
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yJM, 



Figure 9: The ratio Br(uR — » X2) I Br^u^ 



Parameters 


Resolving Methods 


li < M U M 2 


m u ax (xt -> -> Xl) + m T2 


M 2 < Mi 


mu ax (xt -> ^ Xi) 


M 1 < n< M 2 


branching ratio + charge asymmetry 


rrie L <-> m e - fi 


charge asymmetry 



Table 12: Possible flipping solutions and resolving methods. 



mass splitting m-+ — m~+ will be close to m^o — m^o and hence endpoint position 
from (xt ~ * — > Xi ) wm generally differ from m^ ax ~ 112 GeV. 

We end this section by a table summarizing possible flipping solutions and meth- 
ods to resolve it, Table [i~2|. 



5. Implication to LFV Processes 

We begin this subsection by showing the dependence of Br(fi ej) on \i parameter 
for three model points, together with current and future experimental bounds, in 



Figure [L(]. In these plots, we took //, M\ and M 2 of our model points A, A2 and 
A3, and assumed flavor-violating parameters as the model with horizontal symmetry 
described in Section II with x = 0.30. Here, we ignore the left-handed slepton 
mixings since we are interested in the cancellation among diagrams. 

In each plot, there exists a region where the branching ratio becomes very small. 1 
This is due to the cancellation explained in Section II. The positions of cancellation 
point are fairly close for all plots due to the degenerate slepton masses, but only 
point A is in the cancellation region and then gives Br(fi — > ej) below the current 



1 The difference of the minimum value of the branching ratio in each plot is merely the numerical 
artifact. 
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Figure 10: Br(fi — > erf) as a function of /j, for points A, A2 and A3 respectively. We 
assume the MSSM with horizontal symmetry given in Section II. In these plots, x = 0.30 
and a vertical line shows the nominal [/, value for each point. 
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Figure 11: Br(/j — > e7) (red), _Br(/x — > 3e) (blue) and .Br(/x — > e; Ti) (black) as functions 
of // for point A. Here, dashed lines are for the experimental bounds. 




experimental upper bound. Precise determination of fi parameter is then important 
for determination of flavor mixing parameter, especially for points A and A2 which 
have similar collider signature but Br(fi — ► er/) differ by two orders of magnitude. In 
Figure [TT] we show Br(fi —>■ ej), Br (fx — > eee) and \i —>■ e conversion rate in Titanium 
(Br(fi — > e; Ti)) as functions of fi for point A. The fi — > e conversion rate from the 
Br(fi — > &f) also shows strong sensitivity to the SUSY parameters at different fi 
values. Because the cancellation takes place in different region, it is important to 
measure LFV in multiple channels so that we do not miss it accidentally. 

We now estimate the error of the SUSY parameters using information from 
Table [|. First of all, we fixed tan/5 = 10 and required that m^o — m^o, m^o — m^o 
and m^o should not differ from their nominal values greater than 1.5 GeV, 2 GeV 
and 7 GeV respectively. We then looked for the allowed values of Mi,M 2 and fi 
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Parameter 


Central value 


Estimated error 




271.33 


+6.89 
-6.81 


Mi 


122.49 


+7.16 
-7.17 


M 2 


230.89 


+6.57 
-6.54 


/^/Mi 


2.215 


+0.084 
-0.078 


M 2 /Mi 


1.885 


+0.063 
-0.057 



Table 13: Central values and 1-cr constrained intervals of /1, Mi and M2 and their ratios after 
constraints from various cascade decay endpoints are imposed as described in the text. 



Mass difference 


Central value 


Estimated error 


mg R Tfip LR 


1.00 


+0.04 
-0.04 


™e L - m H 


2.00 


+0.48 
-0.49 



Table 14: One-c error estimation of mass splitting between selectron and smuon in GcV when 
m e - H — rrijx R = 1 GeV and mg L — mp, L = 2 GeV are assumed. 



near their nominal values. Furthermore, we calculated m e - fl and using rn™ff\ and 
Tn™n% information and then calculated A% 2 defined in Eq. ( |3.4j ) from other calculated 
endpoints, i.e. m^f 1 ', m^ ax , m™ n , and m™^. The central values of fi,Mi and 

M 2 and their ratios are given with their estimated statistical errors within A% 2 = 1 
(1-cr) region in Table [TB|. The absolute values receive sizable errors while the the 
ratios gives a better sensitivity due to correlation among parameters. Comparing 
Figures |TD| and [II], we can see that the size of error of \x is small enough to resolve 



the position of the cancellation in Br (fx — > ej). 2 

Even if \x is in the region where Br(fi — > ery) is canceled, the direct slepton mass 
measurement at the LHC would provide the information on the non-universality in 
the slepton mass matrices. As had been emphasized earlier, mg R will be measured 
precisely. Then if m e - s and different due to the flavor violating effect, we will 

measure each of them rather accurately. In addition, the error in m e -^ — m^o in Table 
|7] tends to cancel when we take the difference rrie R{L) — rn^ R{L) ■ We therefore estimate 
errors of slepton mass difference in Table by assuming m e - s — = 1 GeV and 
mg L — mp iL = 2 GeV. 

By taking /z, Mi and M 2 of our model point A and assuming flavor-violating 
parameters as the model with horizontal symmetry described in Section II, mass 
difference mg R — Ta^ R and Br(fi — > ej) are plotted as a function of parameter x 
in Figure [12]. Note that since LFV masses in the left-handed slepton sector are 
highly suppressed in this model, then left-handed smuon-selectron mass splitting 
is predicted to be undetectable for this case. The current experimental bound by 
MEGA |l9j and the MEG sensitivity [E2| are also displayed. The mass difference of 



2 Since information of /1 is obtained mainly from m^o , it would not be so sensitive to tan/3. In 



addition, the cancellation point of Br(fi — > &y) , Eq. (2.17) is insensitive to tan/3 neither 
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0.02 



1e-11 
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0.015 



0.01 



- 0.005 



Figure 12: Mass splitting A(m^) 12 /m^ = 2(mg R — m^ fl ) / '(rng R + rnp, R ) and Br([i — ► 
erf) in the model with horizontal symmetry as a function of parameter x. The current 
experimental bound and the MEG sensitivity for \i —* are displayed by horizontal lines. 
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Figure 13: Mass splittings, {rriz L 
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and (mt 



and Br{r — > ^7) in the 



MSSM with right-handed neutrinos as a function of common right-handed neutrino scale. 
h/2L is the lighter/heavier mixed state of left-handed smuon and stau. Dashed lines are 
BABAR and Belle bounds on Br{r — > fj,j). 



order one percent, corresponding to 1 ~ 2 GeV, is allowed by MEGA bound due to 
the partial cancellation. Even though the statistical error of m e - fl — m^ R in Table 14 
are very good, however from the fact that slepton decay width is around 0.2 GeV, we 
simply expect that the LHC has roughly equal sensitivity to the MEG experiment. 



A similar plot for the MSSM with right-handed neutrinos is shown in Figure 13 
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In this model, the left-handed smuon and stau are mixed and form the lighter/heavier 
state h/2L- The horizontal axis now represents common right-handed neutrino mass 
scale. We show Br(r —>■ /ry) because it gives the most stringent bound in this 



model. The experimental bounds of BABAR and Belle [57] allow mass splitting 
between left-handed selectron and In of order few GeV. 

As mentioned before, point A is in the cancellation region and simultaneously 
provides an acceptable dark matter relic density. This can be understood by the 
followings. In Figure [TJ, we plotted values of /x, m^ R , m e - t , m^o, and m^o as functions 
of tuq when min is fixed at 300 GeV and tan/3 = 10. In the plot, we assume the 
mSUGRA relation among slepton and gaugino masses: 



Mi 

2 



0.4mi/ 2 , M 2 = 0.8m 1/2 



m 



m r 



ml + 0.5m 2 / 2 



- — sin 2 9w ) friz cos 2/3, 



(5.1) 
(5.2) 

uiq + 0.15m 2 , 2 — sin 2 9wrn 2 z cos 2/3, (5.3) 

and the value of \x is obtained from the cancellation condition shown in Eq. ( |2.17| ). 
The condition requires \l ~ fnj . Therefore as mo increases, \x increases as well. 
The DM density also increases because Higgsino component of LSP is reduced and 
scalar masses are increased simultaneously. When m reaches value around 150 GeV, 
^dm^ 2 ~ 0.37. In addition, the decay x® ~^ * s no longer open and the slepton 
information may be lost. On the other hand, in a small /i region, DM relic abundance 
is small and Xa &l and x 2 ~^ are always kinematically allowed. In the case, 
the possibility to access to both right- and left-handed slepton masses at the LHC is 
opened if their couplings are not so small. In the same plot, a cross mark represents 
actual \x value for mSUGRA point at mo = 100 GeV. Besides having too large DM 
abundance, it is very far from the cancellation point for Br(fi — > ey). By reducing 
\i from mSUGRA scenario, the acceptable relic density and fi — > ey cancellation 
conditions can be met concurrently. 

Finally we comment on the dark matter 
direct detection cross section of our model 
points. The DM density and spin indepen- 
dent LSP-proton scattering cross section, 
azz, are evaluated for points A, A2 and A3 







a^(10- 8 pb) 


A 


0.1179 


1.55 


A2 


0.0817 


3.15 


A3 


0.0096 


17.50 



Table 15: Dark matter relic density and 
spin independent LSP-proton scattering cross 
section for points A, A2 and A3. 



(Table |15D. Points A2 and A3 are ruled out 
by relic abundance if only one species of DM 
is assumed. Moreover, Point A3 seems to be 
ruled out already by the recent direct detec- 
tion bound from CDMS, a~£ < 4.6 x 10~ 8 pb ||55|| . However, uncertainty in a~£ is 



large due to three important sources. The first one comes from nucleon matrix el- 
ement of strange quark, i.e. (p\m s ss\p), whose recent calculation was found to be 
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Figure 14: Various sparticle masses and \x value as functions of rag at fixed m^ 2 
300 GeV. A cross mark represents fi value of mSUGRA point at mo = 100 GeV. 



an order of magnitude smaller than the previous ones [pq| . Another two sources are 
tan (3 and pseudo-scalar Higgs mass, tua- The pseudo-scalar Higgs exchange diagram 
gives important contribution proportional to j ■ In our analysis, we just take 

tan p = 10 and do not fix tua- However, if the pseudo-scalar Higgs is heavier or tan (3 
has a smaller value than our sample point, the a^ would become smaller and evade 
the direct detection bound. 

We also show in Table |16| the 
sensitivity of Odm ^ 2 and a^ to /i/Mi 
ratio for point A. We assumed NUHM 
with m Hu ^ m and changed 
Table 16: Dark matter relic density and spin inde- value SO that /i/Mi deviates from 
pendent LSP-proton scattering cross section within l-a j ts central value by l-a and 3-cr re- 

and 3-cr deviations of the a/Mi from its nominal value. , mi , • , r 

. spectively. lhe uncertainty of or* 

Here, we assume point A. PX 

is larger than that of Odm« unless 
we fix rriA and tan/3. Although we assume rather good determination of fi/M\ ratio, 
the uncertainty in is not small, i.e. nearly factor 1.5 for l-a. 



/i/M 1 




<J(10- 8 pb) 


l-a 


0.1101 - 0.1271 


1.37 - 1.76 


3-a 


0.0929 - 0.1498 


1.02 - 2.33 



6. Discussion 



In this paper, we studied LHC signature of the one-parameter-extended NUHM, 
7^ n^H D = f»o. The choice of the boundary condition allows the low energy mass 
spectrum Mi < M 2 < A* ~ m i/2- In that case, the LSP relic density may be con- 
sistent with cosmological and astrophysical observations because pair-annihilation 
cross section of the lightest neutralino will be enhanced by the Higgsino compo- 



- 29 - 



nents. We are especially interested in the region where cancellation among leading 
contributions to Br(fi — > cy) occurs in the models with right-handed LFV masses 
because the prediction of Br(fi — > ej) depends strongly on the EW parameters. 
This cancellation occurs in mj ~ /i region, therefore, we take a model point with 
m^o < rrie R < m^o < < m^o as an example. Both the left- and right-handed 
sleptons can be directly produced via neutralino decays at the model point. 

We investigated how well SUSY parameters can be determined at the LHC for 
this choice of parameters. In the region when M\ < M2 < ~ ^1/2 > Xa and X2 
are mixed states with rather large Wino components. They are frequently produced 
from cji, decay and their decay modes into or have large branching ratios. If 
kinematically allowed, would decay dominantly into Wino-like x\\ however, x\ 
and Wino-like X2 a ^ so have small Bino component and they could decay into e#. 
Accordingly, various decay patterns shown in Eq. ( |2.1| ), are expectable at the LHC, 
allowing precise mass determination using endpoint method. 

However, even all above decay could be measured at the LHC and masses of all 
sparticle involving in the decay chains would be identified, there are different regions 
of SUSY parameter space which have mass spectrum consistent with measured end 
points. This ambiguity mostly occurs in the neutralino sector. 

For our model point, we find three solutions with similar mass spectrum but with 
different ordering of /1, Mi and M2 when the relation Mi < M2 is kept. These three 
points predict different Br(fi — > ej) for the same LFV slepton masses. Because the 
original point is so close to the cancellation point, the prediction could differ by factor 
of 0(100) among the three points. The thermal relic density Qdm^ 2 and also 
differ by O(10). However, we find that the reordering leads to different properties of 
neutralinos which will clearly reflect in their decay branching ratio. We showed that 
Br(xi — > 21 + X), rate of events in two hard jets+missing channel, and charge 
asymmetry play important role to lift the degeneracy. This is an excellent example 
of the complementarity of the LHC and the rare decay searches in SUSY studies. 
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A. Appendix 

A.l Trilepton and Four-lepton Distribution in Four-lepton Events 

We have shown in Section III that the four-lepton events from ~> &l ~~ > X2 ~~ > 
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eR — > xi is a useful mode particularly when -Br(%4 — ► ei 



xl) 



is sizable but its 



mu edge can not be seen in two-lepton events. In this Appendix, we show and 
mm' i' distributions as a consistency check. 

Again, we selected four-lepton events by the cuts: 

• exactly two pairs of OSSF leptons where each lepton has p l T > 10 GeV and 
\v\ < 2-5, 

• four leptons must be composed of exactly one e + e~ pair and one pair. 



The cascade decay is 



e L 




Entries 

Mean 

RMS 



xl 



xl 



e R 



e R . The 



The 



1818 
113.3 
46.88 




For each event which passes the cuts, we calculated %f. 
expected to have two trilepton endpoints: 

(1) One from the upper part of the decay chain, i.e. x°a _ 
expected endpoint m^ x = 192.8 GeV. 

(2) Another from the lower part of the decay chain, i.e. ei 
expected endpoint m™fi K = 123.9 GeV. 

By utilizing the m™ ax ~ 123 GeV from xl - * 
ei — > X2 as m Section III, one can further find 
additional constrain on sparticle masses by look- 
ing at the endpoint of the muv distribution as 
follows. Principally, 

(1) if one requires mu < 70.7 GeV and 71 GeV < 
myy < 123 GeV, one expects to see the trilepton 
endpoint from the lower part of the decay chain, 
or 

(2) if one requires myy < 70.7 GeV and 71 GeV < 
mu < 123 GeV instead, the trilepton endpoint 
from the upper part of the decay chain is ex- 
pected to show up. 

However, in our Monte Carlo study, the end- 
point from the upper decay chain is invisible. 
This may be because the distribution for the upper endpoint spreads over wider 
range so that its height of the peak is lower, and then it is buried under back- 
grounds. On the contrary, the one from the lower decay chain in Figure |15| is rather 
impressive as it may receive contribution from chargino decay. 

Finally, mmqi distribution is shown in Figure ITB. The edge is expected not to 
exceed the mass difference m^o — m^o = 208.5. In the left panel, mu of one OSSF 
pair is required to be < 70.7 GeV and that of other OSSF pair is requred to be 
in the range 71 GeV < myy < 123 GeV. Once we further impose the constrain 
may < 124 GeV, the tail almost disappears as shown in the right panel. 



100 150 200 250 

(41 events) 



Figure 15: The trilepton invari- 
ant mass muy distribution (in GeV) 
which satisfies mu < 70.7 GeV and 
71 GeV < myy < 123 GeV. 
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Entries 664 
Mean 179.S 
RMS 31.58 
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Figure 16: The four-lepton invariant mass mui'i> distribution (in GeV) satisfying ran < 
70.7 GeV and 71 GeV < myy < 123 GeV (left panel). The right panel shows the same 
distribution when mm' < 124 GeV is further imposed. 
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